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ABSTRACT 

We report on phase-referenced 23 GHz Very-Long-Baseline-Interferometry (VLBI) observations of the type lib supernova SN 2008ax, 
made with the Very Long Baseline Array (VLBA) on 2 April 2008 (33 days after explosion). These observations resulted in a marginal 
detection of the supernova. The total flux density recovered from our VLBI image is 0.8±0.3 mly (one standard deviation). As it 
appears, the structure may be interpreted as either a core-jet or a double source. However, the supernova structure could be somewhat 
confused with a possible close by noise peak. In such a case, the recovered flux density would decrease to 0.48±0. 12 mly, compatible 
with the flux densities measured with the VLA at epochs close in time to our VLBI observations. The lowest average expansion 
velocities derived from our observations are (1.90 ± 0.30) x 10 s km s _1 (case of a double source) and (5.2 ± 1.3) x 10 4 km s _I (taking 
the weaker source component as a spurious, close by, noise peak, which is the more likely interpretation). These velocities are 7.3 and 
, 2 times higher, respectively, than the maximum ejecta velocity inferred from optical-line observations. 
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j_i 1. Introduction this supernova is qualitatively similar to that of SN 1993J. Since 

supernova SN 2008ax was also cataloged as type lib (Chornock 

Supernova SN2008ax was discovered in galaxy NGC 4490 on 20083, as was SN 1993J, there was evidence to believe that the 

3 March 2008 (Mostardi et al. 2008} at the position a = radio emission of SN2008ax would continue its evolution in a 

12 /! 30"'40.799* and 6 = 4T°3 8'14.8 25". The host galaxy is similar way to SN 1993J's. In such a case, the flux density of 

~8 Mpc distant (de Vaucouleurs l 1 9761 we assume an uncertainty SN 2008ax should have risen well above the VLBI detectability 

of 1 Mpc for the distance), and is dynamically interacting with limit near the end of March 2008. 
NGC 4485 (forming the pair Arp269). Likely as a consequence 
of this interaction, NGC 4490 has a relatively high star forma- 
tion rate (Viallefond et al. 119801 ), which should results in a cor- 
respondingly high supernova rate. 

We assume the discovery date of the supernova was the same On 1 1 March 2008, we proposed a target-of-opportunity 

as the explosion date, since 6 hours before its discovery, the lo- set of global VLBI observations of SN 2008ax at 23 GHz, in or- 

cation of the supernova was imaged with a limiting magnitude of der to detect the supernova radio structure and, possibly, its ex- 

18.3 and no emission was detected (Nakano 2008). Radio emis- pansion. Only antennas of the VLBA were allocated and only 

sion from SN2008ax was monitored with the VLA beginning on 2 April 2008. Unfortunately, the flux density started to drop 

on 3 March 2008, shortly after its discovery (see Stockdale et faster than expected by that time, resulting in only a marginal 

al. [2008a 2008b), with positive detections made at 4.9, 8.4, and detection of the supernova. In the next section we describe the 

22.5 GHz, beginning on 7 March 2008. These detections were all details of our VLBI observations, and in Section [3] we present 

in the millijansky range. The early part of the radio light curve of our results and conclusions. 
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2. Observations and Data Reduction 

We observed supernova 2008ax on 2 April 2008, with the 
VLB A (10 identical 25 m diameter antennas spread over the 
USA from the Virgin Islands to Hawaii). The recording rate 
was set to 256 Mbps, with 2-bit sampling and single polariza- 
tion mode (LCP), covering a total bandwidth of 64 MHz (8 
baseband channels, of 8 MHz width each). Our observations 
were cross-correlated at the Array Operations Center of the 
National Radio Astronomy Observatory (NRAO) in Socorro 
(New Mexico, USA), using an averaging time of 2 seconds. 

The observations of SN2008ax were made in phase- 
reference mode. Each scan of the supernova was of ~2 min- 
utes duration, and short observations (~40 seconds) of strong, 
close by sources were interleaved between these scans of the 
supernova. Since the 23 GHz flux densities of these close by ra- 
dio sources were unknown at the time of the observations, we 
chose the two closest, which were also the strongest at 15 GHz. 
Each pass of the recording tapes (22 minutes long) was then as- 
signed to one of these two calibrators in an alternating scheme. 
The 12-hour long set of observations could, thus, be divided into 
two sets of roughly equal size. In the first one, we observed the 
supernova using the source J1224+4335 as the phase calibrator 
(located 2.23 degrees from the supernova) and in the second one 
we used the source J1225+3914 as the phase calibrator (located 
2.57 degrees from the supernova). 

After the cross correlation, the data were imported into the 
NRAO Astronomical Image Processing System (aips) for cal- 
ibration. We performed the amplitude calibration using gain 
curves and system temperatures measured at all antennas. We 
then used in the scans of the supernova the time-interpolated an- 
tenna gains obtained from hybrid mapping of the calibrators. The 
phase calibration (with account taken of the structures of the cal- 
ibrators) was performed with standard phase-reference calibra- 
tion techniques. The data were then exported for further reduc- 
tion in difmap (Shepherd et al. 1 19951 >. 

3. Results and Conclusions 

The flux density of the calibrator source J1224+4335, obtained 
from hybrid mapping, is 198+1 mJy and the flux density of the 
other calibrator source, J1225+3914, obtained with the same 
procedure, is 122±lmJy. Uncertainties are 3 times the root- 
mean-square, rms, of the hybrid-map residuals (see Readhead 
& Wilkinson 1978). There is no clear detection of SN 2008ax in 
any of the phase-referenced images obtained using both calibra- 
tors. There is not even a clear flux density peak in the images 
obtained from the totality of supernova visibilities (i.e., dynamic 
range above ~6), regardless of the sky coverage of the image or 
the weighting scheme applied in Fourier space. However, when 
the phase-reference calibrator J1224+4335 (i.e., the strongest 
and closest calibrator) is used, there is a possible detection of 
the supernova. The detection arises from applying a visibility 
weighting in Fourier space with the weight of a pixel propor- 
tional to the square of the signal-to-noise ratio (SNR) of the vis- 
ibilities inside that pixel (i.e., we increase the array sensitivity 
in the Fourier inversion). We additionally taper the visibilities 
using a Gaussian, centered at the origin of the uv-plane, with a 
Full Width at Half Maximum (FWHM) of 500 M^. The image 
obtained with such a visibility weighting has a peak flux den- 
sity located at 0.08 mas to the West and 9.2 mas to the South of 
the position used at the correlator, which was taken from X-band 
VLA observations made on 8 March 2009. Therefore, the peak 
flux density detected is located at a = 12 /, 30'"40.79899* and 
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Fig. 1. Hybrid image of source J1224+4335 obtained from our 
observations. The FWHM of the convolving CLEAN beam is 
shown at the bottom-left corner. 



6 = 41°38' 14.8 1580", with an uncertainty of 0.05 mas, which is 
the size of the interferometric beam divided by 2 times the dy- 
namic range of the image (Thomson et al. 1986). Since this po- 
sition is based on a phase-reference to J1224+4335 at 23 GHz, 
we notice that opacity effects in the jet of this source (see Figure 
[TJ could introduce a systematic shift of several mas in the super- 
nova position (the correlation position of J1224+4335 was taken 
from the VLBA Calibrator Survey at 8.4 GHz; see Beasley et al. 
2002). Decreasing the FWHM of the Gaussian taper, or weight- 
ing each pixel with a higher power of the visibility SNR, results 
in a slightly better detection of the supernova but, due to the 
large decrease in resolution, at the price of a detailed detection 
of radio structure. 

After performing a CLEAN deconvolution in the region of 
the flux density peak, we obtain the image shown in Figure [2] 
The rms of the residuals is 0.087 mJy beairT 1 and the peak flux 
density is 0.395 mJy beairT 1 , so the dynamic range of the image 
is 4.7. The total flux density obtained after a deconvolution us- 
ing the CLEAN algorithm is 0.8 mJy. Remarkably, the flux den- 
sity recovered from our VLBI image is a factor ~ 1 .7 larger than 
the flux densities registered with the VLA by Stockdale et al. 
(2008b ) at epochs close in time to the epoch of our observations 
(0.48 ± 0.12mJy on April 1, and 0.45 + O.lOmJy on April 3). 
This large discrepancy indicates that the source structure shown 
in Figure [2] may be a chance superposition of a marginal detec- 
tion (North) and a prominent noise peak (South), as we explain 
below. The integrated flux density of each of the noise peaks of 
the residual image in a 40x40 mas square around the source is 
less than 0.4 mJy (i.e., less than 50% of the integrated flux den- 
sity of the source). We show this wide-field image in Figure[3] 

If the structure shown in Figure [2] is real, does it correspond 
to a core-jet or is it part of a shell? Or, as we suggest above, can 
this structure be a combination of a marginal detection (North) 
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Fig. 2. CLEAN phase-referenced image of SN 2008ax (see text). 
The FWHM of the convolving beam is shown at the bottom-left 
corner. 



with a stronger noise peak (South) than elsewhere in the map? 
Difficult to say. 

Case 1. Partial Shell-like Structure 

For the radio structure part of a shell, we can compare its 50% 
contour level with that of a shell model convolved with the same 
beam. For a shell model with a fractional shell width of 0.3, 
which is the shell width found for SN 1993J by Marcaide et al. 
( 120091 , the outer radius of SN 2008ax would be 1 . 15+0. 15 mas 
(hereafter, all the uncertainties given are equal to the square root 
of the corresponding diagonal element of the covariance matrix, 
with the errors having been first uniformly scaled so that the 
reduced^- 2 is equal to 1). This size translates into an average ex- 
pansion velocity of (4.8 + 0.8) x 10 5 kms -1 , which is superlumi- 
nal. Indeed, we still obtain superluminal expansion velocities if 
we change the fractional shell width to different, unrealistic, val- 
ues such as 0. 1 (a narrower shell width translates into a smaller 
fitted shell size). Hence, it is unlikely that the radio structure is 
part of an expanding shell. 

Case 2. Double Source 

If we instead fit the visibilities to two point sources, one to model 
the brightness peak (at the North) and the other one to model 
the source extension towards the South, we find components 
of 0.40+0.11 and 0.24+0. llmJy, separated by 0.93+0. lOmas. 
This result translates into an average relative velocity between 
components of (3.90 + 0.62) x 10 s km s , which is superluminal. 
If the two components were moving in opposite directions with 
respect to the explosion center, the average expansion velocity of 
the radiostructure would be (1.90 + 0.31) x 10 5 kms _1 , a factor 
~7.3 higher than the maximum ejecta velocity estimated from 
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Fig. 3. Dirty image of SN2008ax with a sky coverage of 
40x40 mas. Due to the limitations of the FFT algorithm, around 
10% of the visibilities (which correspond to the longest base- 
lines) had to be removed in the Fourier inversion for obtaining 
this wide-field image. 



the optical-line emission of this supernova (~ 2.6 x 10 kms , 
Blondin 2008). This velocity is also much higher than the typical 
expansion velocities of the radiostructures of other supernovae 
(~ 1 - 2 x 10 4 kms _1 ). Hence, the two-point source model is 
also unlikely. 

Case 3. Detection with a Close by Noise Peak 

Perhaps, then, we could consider that the structure shown in the 
map is due to a chance (near) superposition of a marginally de- 
tected radio source and a noise peak. In this case, the radio emis- 
sion would not be resolved and its detection would be even more 
marginal. Fitting a shell model (with a fractional width of 0.3) 
to the northern flux density peak results in a source outer diame- 
ter of 0.25 + 0.05 mas and a flux density of 0.48 + 0.12 mJy. We 
notice that this flux density is consistent with the flux densities 
registered by Stockdale et al. 2008b at the same radio frequency 
and at epochs enclosing that of our VLBI observations. The re- 
sulting average expansion velocity is (5.2 ± 1.3) x 10 4 kms _1 , 
a factor of 2 larger than the velocity inferred from optical-line 
emission. This velocity is also ~3 times larger than the expan- 
sion velocities of the other radio supernovae that were observed 
with VLBI (< 2.0 x 10 4 kms _1 ), and would imply (if we as- 
sume a non-decelerated expansion at least until the epoch of our 
observations) a fractional width of the shocked circumstellar re- 
gion of ~0.5. To obtain this estimate, we assume that the optical- 
line emission comes from a region close to the inner edge of the 
shocked ejecta (see Chevalier & Fransson 1994). This fractional 
width is much larger than that found for SN 1993 J and those 
predicted from different models of type II supernovae (Chevalier 
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1982 ). Considering a decelerated supernova expansion would re- 
sult in even larger fractional-shell width estimates. 

General Remarks 

In short, every model used in our analysis results in a supernova 
size much larger than expected from the optical-line velocities 
of this supernova and the expansion velocities found for all su- 
pernovae that could be imaged with VLBI. The lowest average 
expansion velocity compatible with our VLBI data is a factor 
of 2 larger than the velocity inferred from optical-line emission. 
However, despite the apparent significance of our measurements 
(4 cr for the flux density and 5 cr for the size), we have obtained 
only a marginal detection of SN 2008ax with our VLBI observa- 
tions. 
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